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Freezing is one of the most universal phenomena in nature. The abihty to crystalhze is 
common to quite different systems: systems of simphest atoms as well as of large molecules 
of complicated shape, systems of neutral as well as of charged particles, artificial colloidal 
suspensions of polysterene spheres as well as ionic plasma in white dwarfs. 

The freezing transition may be induced either by cooling or by compression, and the 
problem of measuring melting curves is a classic one, having long standing history. To-day, 
due to the diamond anvil technique some melting curves are known with great accuracy and 
up to Megabar pressures. Nevertheless, the detailed microscopic mechanisms even of the 
classical phenomenon of 3D melting are far from being understood. 

There are new problems which gained much interest during last decades connected with 
the problem of non periodic intermediate or ground states. A system which under normal 
conditions freezes into a crystalline phase can be made to freeze into either a quasicrystalline 
or glassy state by a suitable change in the rate of cooling. Further, there are systems in 
nature which do not transform directly from isotropic liquid to crystalline solid but exhibit 
a number of intermediate phases. Among those are surely some 2D systems and liquid 
crystals. 

It is interesting to notice that although the window glass is widely used for a very long 
time, a great number of experimental investigations on glasses has been realized just during 
last decade. Two discoveries seem me the most fascinating among themB. First, the fact 
that the most part of water in the Universe occurs to be in the glass state, being condensed 
from the gaseous phase at very low temperatures. Second, it is just the transition to a 
glassy state and the slowing of diffusion processes that provide the possibility of alimentary 
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product conservation and the life continuation in the desert. 

There is no general theoretical argument that t her mo dynamically stable solid states must 
have a periodic density distribution (seeli). Evidence for the linear stability of aperiodic 
packing to small few-particle displacements can be obtained from a self-consistent phonon 
theoryi. 

In the density- functional theory (DFT) by Singh et al.@ to the glass transition a first- 
order transition at the negative pressure was obtained. In this connection the papers by 
Stishovi'i on phase transitions in expanded matter should be mentioned. Interesting indi- 
cations can be found also in the theory of spin systems. 

The study of disordered spin systems has become a rich and productive science. It 
seemed some time ago that attention completely shifted to neural networks. However, quite 
recently a number of papers has appeared concerning the possilility of obtaining a kind of 
spin glass regime in initially non-random systems. Glasses, of course, have no disorder in 
their Hamiltonian - the randomness is formed in the transition. The deterministic spin 
models with self-induced quenched disorder represent a convenient candidate to describe 
real glasses. Great expectations can be connected with the using of the methods of spin 
glass theory in the theory of real glasses, particularly if combined with DFT of freezing. 

There is another new aspect of the problem in question, which is connected with a new 
concept which appears in the physics of liquid - solid transition just before 1980: bond ori- 
entational ordering. If we consider crystallization as a broken symmetry transition we can 
see that two distinct broken symmetries distinguish crystalline solid from isotropic liquids: 
translational and rotational ones. These two symmetries are not independent, because ro- 
tating one patch of perfect crystal relative to another clearly disrupts not only orientational 
correlations, but translational correlations as well. A relative translation of the two patches, 
on the other hand, decorrelates translational order, but leaves orientational correlations in- 
tact. In such a way one can obtain an exotic fluid or glass state of matter with extended 
correlations in the orientations of locally-defined crystallographic axes, but with short-range 
translational order. This subtle type of order is called bond-orient at ional order (BOO) (a 



"bond" joining near- neighbor atoms) Anisotropic fluids of this kind are a part of recent 
theories of 2D melting! proposed by Halperin and Nelsonll and Young|i§ and based on the 
ideas of Kosterhtz and ThoulessS. This KTHNY theory predicts that the transition may be 
fundamentally different from that observed in ordinary three-dimensional systems. It was 
found that the transition between two-dimensional solid and isotropic liquid can occur via 
two continuous transitions corresponding, respectively, to dissociation of dislocations and 
disclinations. The intermediate phase with BOO was called the hexatic phase. The proper- 
ties of the hexatic phase are similar to those of a nematic liquid crystal, but with a cluster of 
spherical atoms instead of extended molecules. This theory has strong support from exper- 
iments with electrons on helium, computer simulations of 2D electron systems, experiments 
with polystyrene spheres and others. However, a conventional flrst-order transition between 
the two-dimensional solid and isotropic liquid is also a possibility. 

There are some indications that BOO can exist in three-dimensional fluids and glasses, 
too (Steinhardt, Nelson, Toner, Hess, Mitus and Patashinskii, Ryzhov). It should be empha- 
sized, however, that there is a great difference between 2D and 3D cases. In two dimensions 
the symmetry of elementary clusters is crystallographic so that they can simply order in a 
crystal. In three dimensions the energetically preferable symmetry of clusters for standard 
potentials is usually icosahedraJii, incompatible with 3D crystal. In the frame of the flvefold 
symmetry one can pack only 96 (not inflnite number) particles. This means that the cool- 
ing of three-dimensional analog of hexatic phase can not give the icosahedrons packing into 
regular three-dimensional crystal: the clusters must be destroyed before the crystallization 
takes place. It is easier to obtain such a crystal from usual (not supercooled) liquids, where 
there are crystallographic - hexagonal and cubic - clusters. As to supercooled anisotropic 
liquid with icosahedral clusters - it possibly freezes into a glass phase. We shall discuss the 
problem of BOO later in more details. 

Now it is the time to make a remark about 2D crystals. Strictly speaking there is no 2D 
crystals in the thermodynamic limit: they are destroyed by fluctuations. In 1935 Peierls 
has shown that if the temperature is nonzero then the longwavelength phonons destroyed 



long range order in 2D crystals: the meansquare deviations of atoms from their equilib- 
rium positions increase logarithmically with the size of system and the Bragg diffraction 
peacks become smoothed. Later such absence of two-dimensional long-range order was 
proved by Mermin0 with the use of well known Bogoliubov inequalities theorem) for 

correlators^. Now it is clear that one can discuss the quasilongrange translational order in 
2D systems. The appearence of such an order means that the longrange correlations change 
the character of decay from exponential to algebraic. 

In this report we shall discuss some theoretical aspects of the phenomenon of freezing, 
mainly the density-functional theory of the freezing of 3D and 2D systems. The fascinating 
problem of glass formation remains beyond the scope of this report, although there are 
some indications on the applicability of DFT in this case, too. This report has no intention 
of completeness. The choosen list of results and references is influenced by the authors 
own experience in DFT of freezing and in two-dimensional melting as well as by fruitful 
discussions with high-pressure physisists. 

A great deal of information on the phenomenon of freezing can be found in reviews (see, 
for example, the reviews Ref.&Ei and Les Houches sessionsliHEO of 1988 and 1989 for 3D 
systems andillii for the 2D case ) . 

I shall begin with some very simple facts which are not widely known because the crys- 
tallization problems are usually avoided in courses and in standard manuals and textbooks. 

The crystallization transition is strongly first-order in three-dimensions and is marked 
by large discontinuities in entropy, density and order parameters. The correlation functions 
remain short-ranged near the transition. In spite of this the transition is known to have a 
kind of quasiuniversal behaviour. One can mention the following "universalities". 

1. The Lindemann criterion states that the crystal melts when the ratio of the root-mean- 
squared thermal vibration amplitude to the nearest-neighbour spacing attains a certain 
critical value (~ 0.1). The Lindemann criterion has been found to be reasonably correct 
for a number of simple systems. From computer simulations and indirectly, from 
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experiments, the ratio is indeed seen to be about 0.09 for close packed fee structures 
and about 0.12 for open bee structures. This is rather close to the actual melting point 
and gives a quasi- universal criterion of melting. 

2. A comparable one-phase criterion, known as the Hansen- Verlet criterionil, exists on the 
fluid side of the coexistence curve. During a course of extensive computer simulations 
of simple fluids, Hansen and Verlet found that the amplitude of the main peak of the 
structure factor S{kjn) is approximately constant along the crystallization line. There 
is direct experimental evidence that simple fluids as different as sodium and argon, in 
addition to model systems such as a hard-spheres fluid, the LJ fluid, one-component 
plasma (OCP) etc., all freeze when S{km) — 2.9 ± 0.1. 

3. There is a kind of universality in the form of melting curves. Except some exotic curves 
with maxima in p-T coordinates they are monotonic and have no critical end points. 
There are three variants of tricritical points : Liquid-Gas-Solid; Liduidl-Liquidll-Solid; 
Liquid-Solidl-Solidll. The examples of calculations of such curves with tricritical points 
can be found in the recent papers by Bans and coworkers&0 where the role of repulsive 
and attractive parts of pair potential in crystallization is analyzed in detail. 

4. There is a kind of universality in the behaviour of thermodynamic functions along 
melting curves. It seems that the entropy change on melting is constant for a given 
material and that there is the following trendllllll: 

> ln2 when > 0. 

R K 

5. Closely connected with the previous "universality" is a kind of well known "scaling" 
observed in the course of detailed investigation of the freezing of soft spheres systems 
with inverse power potential: 

u{r) = e(cr/r)", 

where e and a measure the strength and the characteristic length of the interaction. 
The reduced excess thermodynamic properties of the soft spheres depend on a single 



variable which is defined as 

7 = {Qa^){e/kBTf^'' = g*T*~^^''. 

Freezing of soft-spheres fiuids have been extensively studied by MC simulation for 
several values of the exponent n, e.g., n = 9, 6 and 4E3'E3. The two extreme 

cases are: the HS system (n = oo), and one component plasma (OCP) {n = 1). The 
computer simulations have revealed some symmetric trends in the melting properties 
of soft spheres when n decreases from oo to ocpili. The relative volume change on 
melting decreases rapidly with n, while the entropy change per particle on melting is 
relatively insensitive to n. Moreover, the more repulsive systems {n > 7) freeze into a 
close-packed fee structure, while the soft repulsions {n < 7) lead to crystallization into 
a bcc phase. The Lindemann and the Hansen- Verlet criterion are exactly obeyed for 
any n. There is also rather general result by WeeksH who has shown for the case of 
arbitrary dimensions that for systems interacting by means of purely repulsive power- 
law potentials, ~ 7--™^ the change in specific volume Av on melting approaches zero 
as m —>■ d, where d is the dimensionality of the system. 

All these facts are pointing towards the possibility of existing an underlying theory. 
But, in fact, there is no such theory so far. It should be emphasized that although the 
freezing is an example of phase transition but different theoretical methods employed for 
the study of other types of phase transitions cannot be employed in this case: 1) Many 
phase transition theories are based on lattice models: many magnetic transitions, gas-liquid 
transition. For the liquid-crystal transition lattice models have not proved very successful 
because the periodicity imposed by the lattice cannot be separated from the spontaneous 
order that ought to arise upon crystallization. 2) The renormalization group (RG) approach 
has not found much use for studying first-order transitions. 3) Computer simulations have 
given a great deal of information, but again have limitations. The system sizes that can be 
simulated are small, and the time scales are short. The question always remains whether the 
time scale of the simulation is sufficient to reach a true equilibrium between two phases. That 
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is why density functional theories (DFT) of freezing occur to be rather useful although they 
have their own limitations, too. DFT is "good" particularly for hard core potentials and for 
high densities because it is a geometrical theory based on packing picture entering the theory 
through direct correlation function (DCF). However the conventional DFT occurs to fail in 
the case of long range potentials and can not describe the melting of some two-dimensional 
systems. 

Now we shall briefly recall the main points of the conventional DFT of freezing, list some 
of the most interesting recent results in DFT and discuss our results on the amelioration 
of this theory. Particularly, we shall reformulate the basic equations of the DFT as to 
describe 2D melting in terms of distribution and correlation functions. The resulting integral 
equations theory can be considered as a variant of DFT theory appropriate to 2D melting 
and, possibly, glass transition. 

The physical idea behind the DFT is the fact that at the freezing transition the corre- 
lation length is only a few atomic spacings. All phenomena at distances greater than the 
correlation length can be treated in a mean-field approximation. Implicit in this approach 
is an assumption according to which a system is either entirely in the liquid or entirely in 
the ordered phase, where no phase coexistence is permitted. Fluctuations are thought to be 
not of great importance in a completely entropy-driven first-order transition. 

The DFT approach to crystallization is based on the theoremil that the Helmholtz free 
energy F[^(r)] of an inhomogeneous system is a unique functional oi the one-particle density 
g{r), which in a crystalline solid is extremely inhomogeneous. The mathematics of DFT 
is the bifurcation theory for the solutions of nonlinear integral equations for one-particle 
distribution function. This idea first appeared in the papers by Kirkwood and MonroeH, 
TyablikovS and Vlasovi^. 

The first papers on the DFT approach to crystallization problem appeared in 1979- 
198lM. 

To describe freezing in three (and two) dimensions we shall follow our papersii'@@@ 
and use a new formalism - that of classical many-particle conditional distribution functions. 



These functions -F(ri|\l/) give the probabihty of finding a particle at ri in the external field 
\[^(ri) The equations for these functions can be obtained from the non-linear integral equation 
for the singlet distribution function F(ri|\I^) in an external field ^^(ri): 

F(ril^) 

P = exp < -/:(W|^rij + + 2^77/ ^fe+Hri, r^+i; 



exp|-/3^(ri) + + ^^y Sfe+i(ri, 



X F{j:2m---F{vk+mdv2...dvk+i}, (1) 

Here z is the activity, p is the mean number density, S'fc+i(ri, rfc+i) is the irreducible 
cluster sum of Mayer functions connecting (at least doubly) k + 1 particles, /3 = l/ksT and 
T is the temperature. 

This equation was derived for canonical ensemble by Arinstein@ on the base of Bogoli- 
ubov functional method^ and then rederived (for the case of small p) by Stillinger and 
BuffS who have used the diagram technique. The simpliest way to derive it in the case of 
grand canonical ensemble was given by RyzhovS. 

If the external field has the form 

*(ri) = E'^(ri-rD 

fc+i 

where $(r) is the interparticle potential, then the function F(ri|\E') is the probability of 
finding a particle at ri, if s particles are at the points r'J, r° 

F(ri|vl/) = F.+i(ri|r?...r°)^^^-^^^"^'"°'-'"°^ 



s\'-li ■■■1 s 

Here -F<i(ri, ...,rs) is the s-particles distribution function. 
In this case the equation (|l|) takes the form 

pF,+i(ri|r?...rO) 



{s k 
k+1 k>l 

xF,+i(ri|r?...r°)...F,+i(r,+i|r?...r°)rfr2...rfrfc+i}. (2) 
The value of z in general case can be obtained from the normalization condition 

i j F,+i(ri|r?...r°)rfri = l. (3) 



If one takes the derivative of (0) relative to ri, one obtains the equihbrium BBGKY 
hierarchyil 

ksTViFs+iir,, r?, r°) + F,+i(ri, r°, r°)Vi " r°) + 

fc+i 

+ p I Vi$(ri - r2)F,+2(ri, ra, r?, .., r°) c/ra = 0, 

along with the explicit expression for Fs+2 as the functional on Fg+i : 

F,+2(ri, T,, r?, r°) = F,(r?, r°)e-^*(^--)F,+i(ri|r?...r°)F,+i(r2|r?...r°) 

pfe-i r dSk+i{ri, ■■■,rfc+i) ^on 



A:>1 

... F,+i(rfc+i|r°...r°)rfr3...rfrfe+i, (4) 

with 

/(ri2) = e-^*(^-) - 1 

This equation gives the exact closure. However it contains infinite series and integrals 
and one has to use some approximations to exploit it. The same can be said about the Eq.(^ 
itself. It is formally closed although the gain is not obvious: the price is the infiniteness of 
series. 

The free energy functional (FEF) of such inhomogeneous system with the density p(r) = 
pFi(r) has the form: 

J^/ksT = J rfrip(ri)[ln(AV(ri) - 1]- 

1 

— Y.k>i (k+iy. / ■ ■ ■ / 'S'fc+i(ri...rfc+i)p(ri) ■ ■ ■p(rfc+i) dri - ■ ■ drk^i. 

or 

J^/ksT = J dr, p(ri)[ln(A^p(ri) - 1] - J^,4p{r)]/kBT. (6) 

The excess free energy J^ex[p{''^)]/kBT is just the generating functional for direct correlation 
functions 

Cn(,ri...r„j - — — — — — . [{ 

5p(ri)---p(r„) 
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If the external field is the potential of atom held fixed at the origin then the functional to 
be extremized depends only on the one particle distribution function g{r) and we can write 
the Taylor expansion for the excess free energy around the liquid in the following form: 

PAF= [ drg{T)ln^-Y,^ [ c(")(ri, r,)Af?(ri)...Af?(rfc)dri...cirfc, (8) 

where 

A^(r) = ^(r) - Qi 

is the local density difference between solid and liquid phase. 

Integral equation for Ag{r) which extremizes AF is formally closed and nonlinear. The 
bifurcation point of the trivial solution of this equation determines temperature and pressure 
when the nontrivial solution appears. This is the point of absolute instability of liquid phase 
against nonconstant density state formation. It can be obtained through exact linearisation 
of the Eq.(|^). We can obtain 

Viln^(ri) = y (ir2C2(ri,r2)V2^i(r2). (9) 
or, for the Fourier-transforms: 

Sg{k) = -PpS{k)U{k), (10) 
S(k) = 1 + p J dre'^^'^[g{r) — 1]. S{k) has the following form in terms of DCF: 

Now one can see that the poles of the structure factor S{k) (the instability points) given by 
the Eq. ([TT|) are the points of mechanical instability of the system. 

In order to obtain the actual point of the phase transition one must use the conditions 
of equality of chemical potentials and pressures for liquid and solid phases and to compare 
the free energy values. 

The full system of equations to be solved in the DFT contains the nonlinear integral 
equation for the function p(r), obtained as the extremum condition for the excess free energy 
and the equilibrium conditions 
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(Ti = T2); Pi = P2; /ii = /is; (12) 

with /ij - the chemical potential and Pi - the pressure of the i-th phase, written in terms of 
the same functions as in (|^). Namely the chemical potential of the inhomogeneous phase is 

= T7 / dr (13) 



V J 6p{r 

The pressure in the solid phase is 

(5"P 



Ap(ri)---Ap(r„)rfri---rfr„. (14) 

p(r)=P; 



5p(ri) • --bpixn 
where P is the functional of p(r) : 

^ = ^{1 ^rp(r)(l-Ci(r;{p(r)}) + /3J-e.(p(r))}. (15) 

To proceed constructively in the frame of the DFT we must choose a concrete form of 
FEF - a kind of closure or truncating - and we must make an ansatz for the average density 
of the crystal. The importance of such ansatz follows from the fact that we are dealing 
with a theory which is equivalent to Gibbs distribution and one has to break symmetry 
following the Bogoliubov concept of quasiaverages@. Now it is necessary to specify the 
crystal symmetry (e.g.lattice type) and to locate the freezing transition for that particular 
lattice type. Other lattice types can then be studied as well. There is always the possibility, 
though, that a more complex lattice that has not been examined will turn out to be more 
stable. 

Let us demonstrate now how the DFT works in the simpliest case of hard sphere system 
truncating the free energy functional in the hypernetted chain (HNC) approximation. The 
one particle distribution function p(r) in the crystalline phase has the form 

p(r) = pi + Ap(r), i J p(r) dr = p,. (16) 

Here p/ and ps - are the averaged density of liquid and solid phase, respectively, and the 
function Ap(r) contains the term having the crystal symmetry. Let us use the following 
actual form for the density change Ap(r) : 
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Ap(r) = Pi Ek V^kc''^'" = Pi'Po + Pi^{r), ^^^^ 

The sum is over reciprocal lattice vectors and the integral is taken over the elementary lattice 
cell A . 

If J labels the sets of the points in reciprocal lattice with equal coefficients (p^ (from 
symmetry condition), then 

Ap(r) = pi^o + PiY: ^jU^), 6(r) = E e^'^'- (18) 
J keJ 

and (pj{r) are independent order parameters. 

In the HNC approximation c„(ri • • -r^) = 0, n > 3, and we obtain from P/ = Pg the 

following equation for the density change Ap(r): 

(1 - piC2{0)) J Ap{r)dr J C2(|ri - rsl) Ap(ri)Ap(r2)ciridr2 = 0, (19) 

or 

(1 - AC2(0))(^oA - 1<pIp1c2{0) - ]^piY.^jC2{kM = 0, (20) 

where 

mj^^j e}(r)cir, (21) 

Uk) = J C2(r)e^^'''-dr. (22) 

The equality of chemical potentials gives the equation 

(1 + v?o) cxp{-piipoC2{0)) = ^ J dr exp jp; ^ C2(/^j)<^jO(r)| • (23) 
Finally, using the equation 

^ J 0{r)0'{r)dr ^ mj5jj>, 

it is easy to obtain 
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^ 1/1// dr^jjv) exp {pi Ej C2(fcj)y?j6(r)} 
"^'l + fo 1/VJ drexp{pij:j52{kj)vAj{r)} ' ^ ^ 

The equations (pOf), (p3|) ( p^) present the closed system of equations for the quantities 
pi,ipo,(pj - the hquid density at the transition, the density change at the crystalhzation and 
the components of the one-particle distribution function, respectively. 

One can also obtain the entropy change in the form 

Even this simple approach gives rather good results. Here in the table we list the result of 
the calculation for HS system in the case of FCC lattice with two reciprocal lattice vector 
values: ipi and v?2- 





Vi 




p* 


As/k^ 


Simulations (MC) 


0.494 


0.103 


11.7 


1.16 


This approach 


0.494 


0.074 


12.6 


1.29 


GELA 


0.495 


0.101 


11.9 


1.15 


SCELA 


0.508 


0.105 


13.3 


1.27 


WDA 


0.480 


0.141 


10.4 


1.41 


MWDA 


0.476 


0.139 


10.1 


1.35 


ELA 


0.520 


0.090 


16.1 


1.36 



It is worth to notice that the interparticle potential (as well as the temperature) does not 
enter the equations. This means, in particular, that they are valid for many-body forces, 
too. The liquid-phase properties enter the theory through the Fourier-transform of the 
DCF, connected with the hquid structure factor S{k) {S{k) = 1 + p J dre^^^[g{r) — 1]) by 
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The structure factor characterizes the relative position of particles in the liquid and so we 
can speak about the DFT approach as a geometrical one and of packing character. 

Let us emphasize that we have truncated the Taylor expansion at second order. Although 
this is a rather crude approximation it is reasonably successful in studying the freezing of a 
number of systems. 

In work done so far using this perturbation theory either c^^\ri, r2, rs) has been ignored 
completely, or only its long-wavelength behaviour has been included, to account for the 
density dependence of the compressibility. 

Although the truncated perturbation approach has significant limitations it has one 
advantage over all approaches proposed to date. This is the fact that it depends only 
on the structure factor of the liquid at a fixed density in the liquid phase region. This 
is experimentally measurable quantity, and this opens up the possibility of studying the 
freezing of rather complex liquids, for which the interaction potentials are not very well 
known, provided that the experimental data on the structure factor exists. 

If the potential is well known for the given real system or if we are dealing with a 
system subject to computer simulations then the direct correlation function to input is 
usually taken from hypernetted chain (HNC) or Percus-Yevic (PY) approximation or from 
the interpolation closure relation of Rogers and Youngil, 



where p{r) = h{r) — c(r), h{r) = g{r) — 1 and /(r) is a "switching function" chosen to be of 
the form 



The parameter k in /(r) is varied until consistency is achieved between the equations of 
state derived from virial and compressibility routes. Another scheme is the modified HNC 



g{r) = exp[-/5M(r)](l + {exp[p(r)/(r)] - l}//(r)). 



(27) 



/(r) = 1 — exp(— Kr) 



(28) 



(MHNC) scheme of Rosenfeld and Ashcroft0 based on the closure 




(29) 
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where the bridge function B(r) is assumed to be a universal function equal to its hard sphere 
(HS) form B{rj,r) calculated for some effective rj] t] is adjusted to yield thermodynamic 
consistency. 

The major limitation of perturbation approaches is that including even third-order term 
is very difficult. The natural question is whether there exists a non-perturbative approach, 
which might include important contributions from the other terms. 

The DFT views the emerging ordered phase as a grossly inhomogeneous liquid with a 
rapidly varying one-particle density p(r), reflecting the lower symmetry of the emerging 
phase. Given such a scenario, it seems unlikely in spite of its success, that the structure 
of the ordered phase can be approximated by the low-order perturbation expansion with 
a uniform system taken as the zeroth or unperturbed system. Motivated by this fact a 
number of workers have attempted to develop approximate but nonperturbative free energy 
functional (FEF). 

This is the idea behind the set of approaches which use an imaginated, unreal liquid at 
an auxiliary density with some kind of variational properties. 

A very popular variant of the nonperturbation approach uses a uniform liquid which has 
different density than pq, so that chemical potential of the solid is no longer equal to that 
of the liquid. The simpliest version of such an "effective-liquid approximation" (ELA) was 
proposed by Bans and coworkers in0 and chooses the reference liquid density in such a way 
that the first reciprocal lattice vector ki of solids matches the first peak in the structure 
factor of the liquid. Such a reference liquid is more dense than the liquid in coexistence with 
crystal at the same temperature, and an extrapolation into the metastable liquid phase is 
required. This is done using HS perturbation theory to calculate the DCF and the structure 
factor over a large range of densities. 

A second set of approaches to constructing an approximate FEF are referred to as 
"weighted density approximations" (WDA). This approach was due to Tarazonali and 
Ashcroft and coworkersil'S. 

The authors have constructed a FEF in such a way that the free-energy density of an 
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inhomogeneous system at a given point is interpreted as that of a homogeneous system, 
but taken at an auxihary density which depends parametrically on the chosen point. The 
effective density is obtained approximately by weighting the physical density over a given 
point. The resultant weighted density approximation (WDA) thus accounts by construction 
for the short-ranged nonlocal effects present in a real, interacting inhomogeneous liquid at 
the given point. 

A simple version of the WDA [known as modified WDA (MWDA)] emerges^ if one 
considers global free energy per particle instead of the local excess free energy per particle. 
In this case the effective density must be position independent. 

Once a FEF has been specified and an ansatz has been made for the crystal density 
(usually it is the form ([T7| ) or a sum over the direct lattice 

N 

gir) = J2cpir-R,) (30) 

i=l 

and assuming a spherically symmetric Gaussian form for 

Lp{r - Ri) = {a/7iY/^exp[-a{r - Ri)^], ) (31) 

one has to calculate as in the simple example considered. 

After 1979 hundreds of calculations were performed of the melting curves of different 
systems. There is a large number of results obtained by the use of DFT which are in 
excellent agreement with real or computer experiments. It should be noticed that according 
to geometrical nature of the theory higher is the pressure, better is the agreement. The 
progress in the DFT as describing the 3D classical melting is connected mainly with the 
papers by Bans, Hansen, Ashcroft and their coworkers. The details of slightly different 
approaches and the lists of results can be found in the reviewsH^ii. 

Technically, the calculations follow the same scheme, mentioned above. The following 
generalizations should be mentioned. 

1. In the case of charged mixture DFT is formulated in terms of number density 

p^(r) = p+{r)+p-{r), 
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charge density 



p^(r) = p+(r) - p (r), 



number-number DCF 




and charge-charge DCF 




-[cSW-c^Ur)]. 



BarratE^ has obtained the phase diagram (using the second-order DFT and the Gaussian 
parametrization for the sohd density) for a mixture of charged HS of the same diameter a 
and opposite charges known as the restricted primitive model (RPM) for ionic hquids. 

2. Freezing of molecular fluids is described in DFT modified to include the angular 
dependence of potentials and DCF of the molecular liquid and to look for angle depending 
solutions for the density. The number of order parameters increases but using of symmetry 
conditions reduces this number. The number of order parameters is determined indirectly by 
the number of harmonics needed to get proper convergence, and this number increases with 
the anisotropy in the shape of molecules. The new phase transition is the orientational phase 
transition. If the anisotropic part of the potential is large enough and the shape of molecules 
differs strongly from the spherical, one has a set of liquid crystal phases. In the opposite case 
one has a small influence of the anisotropic part on the freezing transition and, possibly, an 
orientational phase transition in solid state which may affect or not the crystalline structure. 
All these problems, including a rich phase diagram of liquid and molecular crystal phases 
have been considered in the frame of DFT approach. 

Again in the case of hard-core molecules (hard ellipsoids of revolution (HER), hard 
spherocylindres (HSC) and hard dumbells (HD)) one obtains very good agreement of DFT 
with simulations. In this case the equivalent of HS potential is (hard body potential HB) 



u(ri2, Oi, O2) 



oo,ri2 < L)(fi2, rii, ri2) 



(32) 



0, otherwise. 
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where D{...) is the centre- to-centre distance between the pair of molecules in contact for 
given orientation of the pair. 

3. In order to distinguish between purely geometrical effects and bonding effects charac- 
teristic of real metals Igloi et al@ have subtracted DCF of a HS fluid from the full DCF of 
the metal which they determined from a perturbation theory. We may note that for the HS 
system where the structure is determined solely by the geometrical necessities of a closed 
packing, the structural differences in the free energy (fee and hep) are extremely small. The 
observed differences in metals are due to Ac{k) which measures the bonding effects. 

4. As an example of the application of quantun DFT theory of freezing one can mention 
the paper by Rick et alM. They have used this theory to calculate the phase diagram of 
^He and '^He basing on second-order DFT and taking for the input the correlation functions 
found from the recently developed quantum MC technique. 

All the investigations mentioned above bring to the conclusion that the DFT of freezing 
is a good enough theory in the case of short-range potentials or high pressures, when all 
potentials become more or less hard body potentials. In the case of long-range potentials 
the predictions of DFT theory differ essentially from real and computer experiments. This 
can be demonstrated using as the example the freezing of soft spheres, in which case the fee 
or the bcc phase was found to be a stable structure independendently of n but depending 
on the approximation used. 

Recently in the frame of DFT a number of interesting results were obtained for systems 
widely under investigation now both experimentally and by use of computer simulations. 
Two such systems are particularly "en vogue" now: colloidal suspensions of polydisperse 
spheres and extreme potential systems with isostructural transitions and the absence of 
liquid phase as in Cqq. 

The geometrical nature of the theory defines its success in the case of mixtures of hard 
core particles. Mixtures exhibit phase diagrams which are much richer than those of one 
component systems including crystals substitutionally ordered or disordered. The relative 
stability of these phases depends on the thermodynamic conditions, concentration ratios 
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and ratios of the sizes of atoms. The freezing transition in polydisperse HS system is a 
very interesting problem particularly because there exists now the experimental realization 
of this simple theoretical model - colloidal particles of different size (10-lOOOnm). 

The DFT treatment of Barrat, Bans and Hansen^ was the first to show that, starting 
from monodisperse limit, the freezing transition of the fluid mixture into a constitutionally 
disordered solid changes from a spindle type (at s=0.94) to azeotropic type and (at s=0.92 
/later 0.875) to eutectic type (s = o\la2). Later the DFT approach o was improved inii 
andS. The same results were obtained in simulations by Kranendonk and FrenkelS. 

The crystallization of the polydisperse hard sphere system with continuous distribution 
of particle size (given by gamma- or Gauss distribution) was investigated by use of DFT 
approach in a number of papers. As early as in 1986 inll (see also^) the critical value of 
polydispersity (;^ 5 — 6%) was obtained giving the limit to the possibility of crystalliza- 
tion. Later the simulations Sand the real experiments confirm the existence of such critical 
dispersity (of a little greater value). 

The other interesting experimental results on colloidal suspensions concern the phase 
diagram of "big" spherical colloidal spheres in the liquid of very small polymers. This systems 
was modelled through additional attraction in the hard sphere potential. The range of the 
attraction reflects the polymer size If the polymer size decreases the temperature of the 
liquid-gas transition decreases, too, and becomes lower than the triple point temperature. 
Consequently, the solid phase and only one fluid phase remain on the phase diagramEi^. If 
the range of the attraction is decreasing more an isostructural polymorphic phase transition 
appears in solid phase (in the simulations and in the DFT approach 000Ei^§). The main 
characteristics of phase transitions and the phase diagrams from the DFT calculations are in 
accordance with simulations (iiFil'iiil@) . Analogous results were obtained for hard spheres 
with the Yukawa tail or for adhesive and sticky sphere^^^. 

This problem - the investigation of the attraction role for the crystallization - is a very 
interesting problem, considered as early as by van der Waals. The modern view on the 
problem can be found in HhI. 
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To summarize the first part of the lecture, one can see that the DFT of freezing is simple 
theory which works well enough in the 3D case when one needs not to take into account the 
fluctuations. 

The properties of 2D crystals are quite different from the 3D case. As we have mentioned 
above, in 30-th Peierls has shown that 2D harmonical lattice cannot exist in the thermody- 
namic limit: the meansquare displacement in harmonic crystals has logarithmic singularirity. 
Later the absence of two-dimensional long-range translational order was proved by Mermin 
with the use of well known Bogoliubov inequalities {i/q^ theorem) for the correlators^. 
In fact, in the paper by Merminlll mentioned above (where the absence of crystallographic 
2D order was proved), it was shown that the true long range order - the orientational order 
of bonds - exists in 2D crystal. More precisely, it was shown that the direction of the vector 
between any two neighbouring atoms at finite temperature is the same as at T = 0. 

The bond orientation correlation function < ■d{vi)d{v2) > remains finite at r — > oo. 
Here is the angle between local crystallographic axis and some axis of ideal lattice. In 
continuum approximation it has the form 

^{x, y) = ]^{dxUy - dyU^). (33) 

Anisotropic liquids with bond-orientational order are considered in the phenomenological 
theories of melting, developed by Halperin and Nelson! and Youngjll and based on the 
ideas of Kosterlitz and ThoulessS. This KTHNY theory predicts that the transition may 
be fundamentally different from that observed in ordinary three-dimensional systems. It 
was found that the transition between two-dimensional solid and isotropic liquid can occur 
via two continuous transitions corresponding, respectively, to dissociation of dislocations 
and disclinations. The low-temperature solid phase is characterized by algebraic decay 
of translational order and true long-range bond-orienational order. Dislocations unbind 
at a temperature Tm into a phase with shot-range translational order, but with algebraic 
decay of bond-orientational order. This intermediate phase is called the hexatic phase. The 
properties of the hexatic phase are similar to those of a nematic liquid crystal, except that 
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triangular lattices melt into a phase with persistent sixfold, rather than twofold order. Paired 
disclinations in the hexatic phase ultimately unbind themselves, driving a second transition 
at a higher temperature Ti into an isotropic liquid. 

The base of the KTHNY theory is the mechanism of breaking of quasi-long-range order 
in 2D systems with continuous symmetry developed by Beresinskii0@ and by Kosterlitz 
and ThoulesJl^S. The wel known example of such system is the classical XY model with 
the Hamiltonian: 

^ = -^ E S,S,-^/rfV(V^)^. (34) 

ip is the angle between the vectors Sj and Sj (i, j - the nearest neighbours). 

At low temperature the quasi-long-range order exists characterized by an algebraic cor- 
relation decay < S(r)S(0) >oc r^^. At higher temperature the correlations decay exponen- 
tially < S(r)S(0) >oc exp(— r/^). The symmetry breaking occurs through the appearence 
of free topological defects - vortices: f (Vip) dl = 27iq. 

The transition temperature can be obtained simply from energetic ballance using the 
fact that the vortex energy is (from (0)) 

J fL 27r 



Ey = — I — dr = j7rln(L/a), 
2 JO r 



where a is the lattice costant and L - the system size. The creation of a vortex changes the 
free energy: AF = Ey~TS^ where the vortex entropy S = 2kB ln(L/a). At the temperature 
T > Tq = TrJ/ks the value of AF = (Jvr— 2A;bT) ln(L/a) becomes negative so that the vortex 
creation becomes energetically profitable. The Hamiltonian for vortex-vortex interaction is 
equivalent to that of 2D Coulomb gas: 



|r — r'' 



H, = -nJ g(r)g(rOln^ + E,Y.<l\^)- (35) 

|r-r'|>a ^ r 

The physics of KT transitionS is the dissociation of vortex pairs in the presence of 
the screening. The KT theory presents a renormgroup approach to screening. The KT 
transition is a continuous transition from low temperature quasi-long-range ordered phase 
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to high temperature disordered phase. The Tkt is just that obtained above from simple 
thermodynamical consideration. 

Based on the KT theory the KTHNY theory describes the 2D melting as two continuous 
KT transitions. The first transition - the dislocation pairs unbinding, the second one - the 
unbinding of disclination pairs. 

The description of the first transition in KTHNY theory is based on the elastic Hamil- 
tonian for 2D triangle lattice: 



2^uf, + Am 



kk 



(36) 



where 



Ui 



dui(r) ^ du^(r) 



(37) 



and /i and A are the Lame coefficients. 

Free dislocations (analogs of the vortices in XY model) break the quasi-long-ranged 
periodic translational order and cause the shear modulus fi to vanish. 

The dislocation Hamiltonian has the form: 



dis 



E {b(r.)b(r,) In - (b(r.)r.,0(b(r.)r.,) | ^ 

M 



+ E,^b2(r,), 

where Ed is the core energy of dislocation. 

^ 4;u(/i + A) 
2/x + A 

The KT phase transition to the hexatic phase takes place at the temperature 



(38) 



(39) 



167r. 



(40) 



Following Halperin and Nelsoni, the phenomenological orientational order parameter is 



^(r) = e^^''^'-), 
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(41) 



where 'd{r) is the bond orientation (|3^). The quasi-long-ranged orientational order at T > 
is described in terms of algebraic correlation decay: 

< ^/'*(r)^(0) >oc r-'^«(^\ (42) 

The Hamiltonian for the hexatic phase has the formH: 

HA = lKAiT) I dV(V^(r))^ (43) 

The Frank constant K^iT) is connected to 1]^ through the equation: 

In the solid phase the disclinations are tightly bound into pairs (the dislocations), how- 
ever in the hexatic phase the interaction is screened by dislocations and so occurs to be 
logarithmici'0, and so for the disclination Hamiltonian in hexatic phase we obtain the 2D 
Coulomb gas like form againHS: 

H,..c = -^^^^ E ^(r)^(r') In + E,, Yl ^'(r), (45) 

where Ed is the core energy of disclination, s(r) = 1 for the atom with 7 nearest neighbours, 
s(r) = —1 for 5 nearest neighbours. 

The second KT transition - the disclination pairs unbinding - takes place at the 
temperature! Tj > T^: 

= (46) 

The KTHNY picture is the following. At T = there exist long-range translational and 
orientational orders, however, at T > the long-range translational order is transformed into 
the quasi-long-range one due to smooth phase fluctuations, the orientational order remaining 
unchanged. At the melting temperature Tm, the quasi-long-range order is destroyed by the 
singular phase fluctuations (dislocations). The appearance of free dislocations means that 
the system ceases to offer any resistance to shear (/i = 0), i.e., it becomes a liquid. It should 
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be emphasized that the amphtude of the order parameter pc does not become equal to 
zero at T^, but it does at the some mean-field temperature Tmf, which can be determined 
by equating the free energies of the solid and liquid phases as the functionals of the local 
density. (It is important to note that it is this point Tmf that has been considered by 
Ramakrishnanii as the melting point.) The second transition into an isotropic liquid takes 
place at a higher temperature Tj . 

This theory has strong support from experiments with electrons on heliumli2^'E23, com- 
puter simulations of the 2D electron systemsEslllii^elasticity simulations^, experiments with 
polystyrene spheresEsl (however, the topological defects in this case are complex and are not 
consistent with a simple KTHNY picture of melting). 

A conventional first-order transition between the two-dimensional solid and isotropic 
liquid is also a possibility. Several theories predict a single first-order melting transition in 
The sum of the experimental evidence on adsorbed atoms is weighted 



against an interpretation in terms of KTHNY melting and seems to show a weak first- 
order transition0. Most simulations studies of 2D melting indicate that the strictly two- 
dimensional hard-core potential systems melt via first-order transition (see, for example 



and references therein), as do studies of systems interacting with intermediate strength 
potentialjnH. 

Therefore, we expect that the melting behavior may depend crucially on the interaction 
potential, and the first-order character of the transition is weakened as the potential is 
softened. 

The simple physical picture may be painted in order to illustrate this possibilityEl 



lie 



As mentioned above, there is an analogy between the hexatic phase and a nematic liquid 
crystal, but the role of the rodlike molecules is played by the hexagon clusters consisting of 
an atom with its nearest neighbors. Ordering of these hexagons is possible only if the range 
of the interparticle interaction is large enough to provide the interaction between clusters. 
Therefore, the interparticle interaction should extend at least over several interparticle dis- 
tances, and hexatic phase could exist for long-range potentials, but not for shot-range ones. 
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However, only first-principles study could put these heuristic reasonings on firmer footing. 

The microscopic theory of 2D melting can be obtained as a generalization of the 
DFT theory of freezing described above. We developed an approach appropriable to 2D 



meltinglii31^li22rilHa'tl3'E13. Our approach can predict, basing on the knowledge of interpar- 
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tide potential, which scenario is to be realized: 1) < Tmf, the system melts by means 
of two continuous transitions of KT type; 2) Tmf < T^, the system melts by means of a 
first-order transition. 

Our approach differs from the standard DFT theory of freezing in two main points: 
First, we permit the Fourier coefficients PG(r) of the one-particle distribution function 
expanded in a Fourier series in reciprocal-lattice vectors {G}: 



P(r) = ^PG(r)e 

G 



Gr 



to fluctuate - to vary slowly over distances of order G ^ and to have the amplitude and 
phase: 

PG(r)=PGe^''"^'-). 

Here u(r) is the displacement field. 

Second, we permit the liquid to be anisotropic: we consider as possible the existence 
of a phase with constant density but angular dependent two-particle distribution function 
i^2(ri - ro) ^ fi'(rio). 

These two points of generalization define the two new order parameters: the fluctuating 
Pcir) and the Fourier coefficients characteristic for the broken symmetry of the function 
F2(ri — ro). Our approach again is based on the Eq.(Q) of previous section. 

Let us describe in details the microscopic characteristics of a hexatic phase. 

The relative spatial distribution of pairs of particles is characterized by the function 
F2(ri|ro) = F2{ri — ro). The vector ri — ro defines the direction of the bond between the 
molecules at the points ri and tq. In the ordinary isotropic liquid the nearest neighbouring 
of a given molecule (the first coordination sphere) has a definite local symmetry, which 
can be characterized by the set of bond directions. The local structure of the liquid in 
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the neighbourhood of a molecule at the point Tq is characterized by the bond directions 
r' = r2 — Tq. It occurs that if the point Tq is at sufficiently large distance from fq then there 
is no correlation between the directions r = ri — and r' = r2 — Tq. In this case after the 
averaging over the system as a whole the pair distribution function transforms into the RDF 
and the equation (0) for s = 1 has the solution F2(ri — tq) = g{\ri — ro\), which corresponds 
to ordinary isotropic liquid. 

When we approach the anisotropic liquid phase the long-ranged correlations between 
the bond directions r and r' do appear and the averaged two-particle distribution function 
depends on the bond direction now. 

In the vicinity of the transition one can write 

F2{r^, To) = g{\r^ - ro|)(l + f{r, - ro)) (47) 

where /(ri|ro) has the symmetry of the local neighbourhood of the particle at Tq. In the 
two-dimensional case 

/(ri|ro) = /(ao,ro,¥5). (48) 

Here o = |ri — rol; c^o is the radius of the first coordination sphere, and (p is the angle of the 
vector ao. The function / may be expanded in a Fourier series 

oo 

/(ao,ro,<^)= /m(ao,ro)e*'"'^. (49) 

m,=— oo 

The Fourier coefficients define the order parameters. These parameters become nonzero 
at the temperature mf defined by the bifurcation condition, that is by the eigenvalue 
of the linearized (relative to /) eq. (0) with s = 1. 

At the same time, when one approaches the line defined by the bifurcation condition, the 
correlation radius for the orientation fiuctuations of the pair distribution function diverges. 
This fact can be shown with the use of the gradient expansion technique in the case of the 
equation (0) for s = 3, if we write the long range part of the correlator using the principle 
of a weakening of correlations (§) as: 
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F4(ri, r4) = g{\vi - rsD^drg - r4|)(l + /4(ri, r4)) 



(50) 



/4(ri, r4) = fi{r, R, p, ^i, ^2)- 

Here is the angle between the vector r = ri — r2 and the axis R = r2 — r3, (y92 is the angle 
between the vector p = — and the same axis. We have f4^{r, R, p,ipi,ip2) — > when 
i? — > 00 . 

The function /4 can be presented as a Fourier series 

/4(r, R, p, (fi, (P2) = fpqi^^ P) exp(zpv5i + iq(p2). (51) 
The asymptotics of the solution is of the form 

Here the correlation length is characterized by the properties of the isotropic liquid and 
diverges in the vicinity of the bifurcation line. The asymptotic behaviour of the correlation 
function /g^e , derived in our approach, (far away from the transition) can be compared with 
that of the phenomenological KT theory, thus giving the microscopic expression for Frank 



Ka{T) = 648fcBTa^|/6|'[r6(ao, ao) - l/2(r5(ao, ao) + rr(ao, ao))], (52) 
1 /"^i 



constan^'^^ 



Zn Jo 



where = 2/(-\/3p). The function r(ri, r2, ^p) was introduced in our paperstl3ilHy on bond 
orientational order and has the form 

r(ri,r2,¥?) = ^ — — / 5fc+i(ri, ...,rfc+i) 

fc>i \^ ~ 

X 5((|r3 - rol) ■ ■ ■ g{\rk+i - roDciVa ■ ■ ■ dVfe+i, (53) 

where, as earlier, 5*^+1 (ri, r^+i) is the irreducible cluster sum of Mayer functions connect- 
ing (at least doubly) k+1 particles, ri = |ri — ro|, r2 = |r2 — ro|, and ip is the angle between 
the vectors ri and V2- 
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In the spirit of DFT we expand the function (|53|) in a functional Taylor series in powers 
of h{r) = g[r) — 1: 

Coo n I- 
n=l 

X /i(|r3 - rol) ■ ■ ■ /i(|r„,+2 - roD^Vg ■ ■ ■ rfV^+a) • (54) 

In the hypernetted chain approximation we can drop in the expansion (^4]) all terms with 
n > 3 and obtain 

r(ri,r2,(^) = pc(2)(|ri-r2|). (55) 
In this approximation, Eq. (|5^) has the form 

Ka{T) = 648/cBTao|/6p[c6(ao, Oq) - l/2(c5(ao, Oq) + 07(00, ao))], (56) 

where 

1 f^TT 



Cm[ao,ao) = — / rfv3c^^)(y2ao(l - cosv?)^''^;p)e" 

/TT JO 



■imip 



Now the equation for the absolute value of the order parameter ie 
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120. 



^ /o^'' dip cos6v9 exp[12c6(ao,ao) cos(6v9)/6] ^^^^ 
Sq"" dip exp[12c6(ao,ao) cos(6v3)/6] 

(We have made use of the Eq.(0) for the function F2(ri|ro)). As in the case of standard 
DFT this equation corresponds to the minimum of the free energy functional. 

The phase transition from the hexatic phase to the 2D solid phase can be analyzed in 
a similar way. In this case the correlator which diverges is the density-density correla- 
tion function, the density Fourier components being the order parameters (pG(r)pG(O)) oc 
^g(i"); h2{r) = J2g hG{r)e''^^, where G are the reciprocal lattice vectors. 

p(r) = pF, (r) = po + Ap(r) = p, + Y^p^ (r) e^^"" (58) 

G 

PG(r) = IpGle*"^'-) 



28 



In both cases the solution with broken symmetry appears at the point where the corre- 
lation decay changes its character. 

The microscopic expressions for the nonrenormalized Lame coefficients were obtained 
ii^M by comparing the long-range behaviour of our correlator with the asymptotic behavior 
of the order parameter correlator from the phenomenological elastic energyi. For short- 
range potentials, expressions were obtained inliHl for the elastic moduli which correspond to 
fluctuations with wave vector equal to the smallest reciprocal lattice vector. Generalization 
to the case of an arbitrary number of inverse lattice vectors leads to the following expression 
for the Lame coefficients: 

/^=^Epg^gG^(7c + 25g), (59) 
16p ^ 

>^ = ^Y.Pl^GG'{jc-65a), (60) 

where 

7G = 271 p I r^drc^^\r-p)Jo{Gr), (61) 
6g = 2'kpJ r^drc^'^\r;p)Ji{Gr)/{Gr), (62) 

Jo{x) and Ji{x) are the Bessel functions, and rriG is the number of reciprocal lattice vectors 
with the same length. 

In the case of the long-range Coulomb interaction, an additional term arises in the elastic 
HamiltonianllHi, which makes the effective modulus A diverge, A = oo, but the expression 
(^)for fi remains to be valid. The modulus K takes the form: 

K = Afial/keT. (63) 



Inlii3 the authors present the first-principles estimates for the stability limits of the solid 



and hexatic phases for the 2D electron system and the system of hard disks which are two 



opposite cases of the 1/r" potential. InUi^ 2D system of vortices in a superconducting film 
interacting via a potential which is even softer than in 2D electron system was considered. 
The transition temperatures and Tj were obtained from the KT theory using microscopic 
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expressions for elastic and Frank moduli and were compared with the value of Tmf from the 
standard DFT theory. 



1. Beginning with the classic paper by Alder and Wainwrightli^ it has been assumed that 



a system of hard disks melts by means of the first order transition. However, only recently 



have convincing proofs of this fact been obtainedlii3. The standard DFT approach gives good 
qualitative and quantitative agreement with the results of computer simulations. In the 
simpliest version we obtain ps = 0.933 (to be compared with 0.921 from simulationsliHi) . At 
the transition point we have K/IGtt = 6.29 which is much larger than the value K/16tt = 1, 
at which the dissociation of dislocation pairs takes place. 

However, an even more convincing argument in favor of the first order phase transition 
is provided by the analysis of the possible existence of hexatic phase. Using the equation 
(p6|) we obtain for the Frank constant for the system of hard disks Ka{T) < at all density, 
so that the hexatic phase cannot exist. 

2. In the case of 2D classical Wigner crystal we have obtained the two-stage scenario 
of melting. In terms of the dimensionless parameter F = (npY^'^e^ / (ksT) our results are: 
Tmf = 21.58, F^^ = 80.15 and Tf^ = 24.5. However, the renormalization of parameters 
depends strongly on the unknown values of the core energies of disclinations and dislocations 
and this can change the result. 

3. In thin superconducting film the interaction energy of two vortices located at the 
points r j and r ^ (r jj = | r j — r ^ | ^ ^ has the formE^i 



Ho (^) - Yo 



$(r„)^ -Jrlnf^l r„«A (64) 



tjJ ~ 47r2A A / 

Hr,,) ^ ^ r,, » A 

where A(T) = 2\\(T)/d is the effective penetration depth, d is the film thickness, is 
the bulk penetration depth, ip^ = hc/2e is the flux quantum, Hq{x) - the Struve function 



and Yo{x) - the Neumann function. The potential (|6^) is long-ranged (it is even "softer" 



than the Coulomb potential 1/r) and one might therefore expect that vortex lattice will 
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melt through two continuous transitions. As to real and computer experiments, not only 
the type of the transition, but even the very existence of vortex lattice melting, remains an 
open question. This system was investigated in details by the authors inE 



1.3C-134 



On the basis of our approach we have made the calculations for the specific case of a 



niobium film of thickness 20A, investigated experimentally in Ref.li^. The region of the 
hexatic phase is very well pronounced on the H-T phase diagram. If the film thickness 
increases this region narrows (see figures inlli^). 
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16211). 



^ See a number of the articles in Science, 267, 1924-1951 (1995). 
^D.Ruelle, Physica A 113, 619 (1982). 

3 J.P.Stoessel, P.G.Wolynes, J. Chem. Phys. 80, 4502 (1984). 

^ Y.Singh, J.P.Stoessel, P.G.Wolynes, Phys. Rev. Lett. 54 1059 (1985). 

^S.M.Stishov, JETP Lett. 57, 196 (1993). 

6 S.M.Stishov, A.V.Smirnov, JETP Lett. 57, 727 (1993). 

"^M.Widom, K.J.Strandburg, R.H.Swendsen, Phys. Rev. Lett. 58, 706 (1987). 

^ For a recent review, see K. J. Strandburg, Rev. Mod. Phys. 60, 161 (1988), and references 
therein. 

^B. I. Halperin and D. R. Nelson, Phys. Rev. Lett. 41, 121 (1978); D. R. Nelson and B. I. 
Halperin, Phys. Rev. B 19, 2457 (1979). 

OA. P. Young, Phys. Rev. B 19, 1855 (1979). 

31 



^1 M. Kosterlitz and D. J. Thouless, J. Phys. C 6, 1181 (1973). 
12 F. C. Prank, Proc. R. Soc. London Ser.A 215, 43 (1952). 
1^ R. E. Peierls, Ann. Inst. Henri Poincare, 5, 177 (1935). 
i^N. D. Mermin, Phys. Rev. 176, 250 (1968). 

15 N.N.Bogoliubov, JINR Preprint R-1451, 1963. 

16 S.M.Stishov, Soviet Phys. Uspechi (UFN) 114, 3 (1974). 
i'^ S.M.Stishov, Soviet Phys. Uspechi (UFN) 154, 93 (1988). 
1^ Ya. Singh, Phys. Rep. (Phys. Lett.C) 207 N6 351, (1991). 

1^ "Simple molecular systems at very high density", Les Houches (1988) A.Polian, P.Loubeyre and 
N.Boccara editors, Plenum Press, N.Y. and London. 

2^ "Liquid, freezing and glass transition", Les Houches (1989) J. P. Hansen, D. Levesque, J. Zinn- 
Justin editors. North Holland. 

21 For a recent review, see K. J. Strandburg, Rev. Mod. Phys. 60, 161 (1988), and references 
therein. 

22 K.J. Strandburg (editor), Bond-Orientational Order in Condensed Matter Systems. (Springer- 
Verlag, New York, 1992). 

23 J.RHansen, L.Verlet, Phys. Rev. 184, 150 (1969). 

2^0. F. Tejero, A.Daanoun, H.N.W.Lekkerkerker, M.Baus, Phys. Rev. Lett. 73, 752 (1994). 

25 A.Daanoun, H.N.W.Lekkerkerker, M.Baus, Phys. Rev. E 50, 2913 (1994). 

26 c. F. Tejero, A.Daanoun, H.N.W.Lekkerkerker, M.Baus, Phys. Rev. E 51, 558 (1995). 
2^^ J. P. Hansen, Phys. Rev. A 2, 221 (1970). 

2^ W. G. Hoover, M. Rose, K. W. Jonson, D. Henderson, J. A. Barker, B. C. Brown, J. Chem. 

32 



Phys. 52, 4931 (1970). 

29 W. G. Hoover, S. G. Gray, K. W. Jonson, J. Chem. Phys. 55, 1128 (1971). 

30 J. P. Hansen, D. Schiff, Mol. Phys. 25, 1281 (1971). 

31 W. G. Hoover, M. Ross, Contemp. Phys. 12, 339 (1971). 

32 W. G. Hoover, D. A. Young, R. Grover, J. Chem. Phys. 56, 2207 (1972). 

33 J.D.Weeks, Phys. Rev. B24, 1530 (1981). 

34 R. Evans, Adv. Phys. 28, 143 (1979). 

35 J.G.Kirkwood, E.Monroe, J. Chem. Phys., 9, 514 (1941). 

36 S.V.Tyabhkov, JETP 17, 386 (1947). 

3''' A.A.Vlasov, J. of Phys. 9, 130 (1945); Many body system theory, Moscow, Gostehizdat, 1950. 

38 T.V.Ramakrishnan, M.Youssouff, Phys. Rev. B 19, 2775 (1979). 

39 V.N.Ryzhov, E.E.Tareyeva, Phys.Lett.A, 75, 88 (1979). 

40 V.N.Ryzhov, E.E.Tareyeva, Theor.Math.Phys.(Moscou), 48, 416 (1981). 

41 A. D. J. Haymet and D. W. Oxtoby, J. Chem. Phys. 74, 2559 (1981). 

42 V.N.Ryzhov, E.E.Tareyeva, Physica A, 109A, 357 (1981). 

43 V.N.Ryzhov, Theor.Math.Phys. (Moscow), 55, 128 (1983). 

44 E.A.Arinshtain, Dokl. Akad. Nauk (USSR), 112, 615 (1957). 

45 N.N.Bogohubov, Problems of dynamical theory in statistical physics (Moscou, Gostehisdat, 
1946). 

46 F.H.Stilhnger, F.P.BufF, J.Chem.Phys., 37, 1 (1962). 
4'^ V.N.Ryzhov, thesis (JINR, Dubna, Russia, 1989). 

33 



N.N.Bogoliubov, JINR Preprint R-1451, Dubna, 1963; Phys.Abh.S.U., 6, 1, 113, 229 (1962). 
''^F. J. Roger, D. A. Young, Phys. Rev. A 30, 999 (1984). 
^°Y. Rosenfeld, N. W. Ashcroft, Phys. Rev. A 20, 1208 (1979). 
s^M.Baus, J.L.Colot, Mol.Phys.55, 653 (1985). 

52 P.Tarazona, Mol.Phys. 52, 81 (1984); Phys. Rev. A 31, 2672 (1985). 

53 W.A.Curtin, N.W.Ashcroft, Phys. Rev. A32, 2909 (1985). 

54 A.R.Denton, N.W.Ashcroft, Phys. Rev. A39, 4701 (1989). 

55 A.R.Denton, N.W.Ashcroft, Phys. Rev. A 39, 426 (1989). 

56 Y. Singh, Phys. Rep. 207, 351 (1991). 

5^M. Bans, J. Phys.: Condens. Matter 1, 3131 (1989). 

58 J.P.Hansen inS 

59 D.W.Oxtoby in§ 

6°H.L6wen, Phys.Rep. 237, 249 (1994). 

61 J.L.Barat, J. Phys. C 20, 1031 (1987). 

62 F. Igloi, G. Kahl, J.Hafner, J. Phys. C 20, 1803 (1987). 
63S.W.Rick, J.D.McCoy, A.D.J.Haymet, J. Chem. Phys. 92, 3040 (1990). 

64 J.L.Barrat, M.Baus, J.P.Hansen, Phys. Rev. Lett. 56, 1063 (1986); J. Phys. C 20, 1413 (1987). 

65x.C.Zeng, D.W.Oxtoby, Phys. Rev. A 41, 7094 (1990). 

66 A.R.Denton, N.W.Ashcroft, Phys. Rev. A 42, 7312 (1990). 

6^ W.G.T.Kranendonk, D.Frenkel, Mol. Phys. 72, 679, 715 (1991). 

68 J.L.Barrat and J.P.Hansen, J.Phys.46, 1547 (1986). 

34 



69R.McRae, A.D.J.Haymet, J. Chem. Phys. 88, 1114 (1988). 
^° P.G.Bolhnis and D.A.Kofke, Phys.Rev.E 54, 634 (1996). 

R.T.Farouki and S.Hamaguchi, J.Chem.Phys. 101, 9885 (1994). 
^^E.L.Pollok and J.P.Hansen, Phys.Rev. A 8, 3110 (1973). 

J.Medeiros e Silva and B.J.Mokross, Phys.Rev. B 21, 1972 (1980). 
^^H.Totsuji and K.Tokami, Phys.Rev. A 30, 3175 (1984). 

75 W.L.Slattery, G.Doolen and H.E.DeWitt, Phys.Rev. A 21, 2087 (1980); ibid 26, 2255 (1982); 
G.S.Stringfellow, H.E.DeWitt and W.L.Slattery, ibid 41, 1105 (1990). 

76 M.O.Robbins, K.Kremer and G.S.Grest, J.Chem.Phys. 88, 3286 (1988); K.Kremer, 
M.O.Robbins and G.S.Grest, Phys.Rev.Lett. 57, 2694 (1986). 

77 C.Rascon, L.Mederos and G.Navascues, J.Chem.Phys. 103, 9795 (1995); Phys.Rev. B 51, 14899 
(1995). 

7SE.J.Meijer and D.Prenkel, J.Chem.Phys. 94, 2269 (1991). 

79 M.J.Stevens and M.O.Robbins, J.Chem.Phys. 98, 2319 (1993). 

s°G.Dupont, S.Moulinasse, J.P.Ryckaert and M.Baus, Mol.Phys. 79, 453 (1993). 

M.H.J.Hagen, E.J.Meijer, G.c.A.M.Mooji, D.Prenkel and H.N.W.Lekkerkerker, Nature (Lon- 
don) 365, 425 (1993). 

82 C.F.Tejero, A.Daanoum, H.N.W.Lekkerkerker and M.Baus, Phys.Rev.lett. 73, 752 (1994). 

83 L. Mederos and G. Navascues, Phys. Rev. B 50, 1301 (1994). 
8^T. L Schelkacheva, Phys. Lett. A 214, 95 (1996). 

85 M.Hasegawa, K.Ohno, Phys.Rev. E 54, 3928 (1996). 

86 C.Rascon, L.Mederos and G.Navascues, Phys.Rev. E 53, 5698 (1996); Phys.Rev. E 54, 1261 

35 



(1996); 

^^B.B. Laird and D.M. Kroll, Phys. Rev. A 42, 4810 (1990). 

88 A. Kloczkowski and A. Samborski, J. Chem. Phys. 88, 5834 (1988). 

89 S. J. Smithline and A. D. J. Haymet, J. Chem. Phys., 83, 4103 (1985). 

90 C.F. Tejero and M. Bans, Phys. Rev. E 48, 3793 (1993). 

91 B.Bagchi and C. Cerjan, Phys. Rev. A 31, 1647 (1985). 

92 J.D.Weeks, R.L.B.Sehnger and J.Q.Broughton, Phys.Rev.Lett. 75, 2694 (1995). 

93 ... 1996 . 

9^ . . , 59, 907 (1970). 

95 . . , 61, 1144 (1971). 

96 J. M. Kosterhtz, J. Phys. C 7, 1046 (1974). 

97 D. R. Nelson and R. A. Pelcovits, Phys. Rev. B 16, 2191 (1977). 

98 D. R. Nelson, Phys. Rev. B 26, 269 (1982). 

99 T. V. Ramakrishnan, Phys. Rev. Lett. 48, 541 (1982). 

C. C. Grimes and G. Adams, Phys. Rev. Lett. 42, 795 (1979). 

101 R. Mehrotra, B. M. Guenin, and A. J. Dahm, ibid. 48, 641 (1982); F. Gallet, G. Deville, A. 
Valdes, and F. L Williams, ibid. 49, 212 (1982); G. Deville, A. Valdes, E. Y. Andrei, and F. 
L WilUams, ibid. 53, 588 (1984); R. Mehrotra, C. J. Guo, Y. Z. Ruan, D. B. Mast, and A. 
J. Dahm, Phys. Rev. B 29, 5239 (1984); M. A. Stan and A. J. Dahm, Phys. Rev. B 40, 8995 
(1989). 

102 R. H. Morf, Phys. Rev. Lett. 43, 931 (1979). 

103 V. M. Bedanov, G. V. Gadiak, and Yu. E. Lozovik, Zh. Eksp. Teor. Fiz. 88, 1622 (1985) [Sov. 

36 



Phys. JETP 61, 967 (1985)]. 
i°^Y. Saito, Phys. Rev. B 26, 6239 (1982). 

c. A. Murray and D. H. Van Winkle, Phys. Rev. Lett. 58, 1200 (1987). 
s. T. Chui, Phys. Rev. B 28, 178 (1983). 

107 H. Kleinert, Phys. Lett. 95A, 381 (1983); 

108 V. N. Ryzhov, Theor. Math. Phys. (Moscow) 88, 449 (1991). 

109 V. N. Ryzhov, Zh. Eksp. Teor. Fiz. 100, 1627 (1991) [Sov. Phys. JETP 73, 899 (1991)]. 

110 Yu. E. Lozovik and V. M. Farztdinov, Sohd State Commun. 54, 725 (1985); Yu. E. Lozovik, V. 
M. Farztdinov, B. Abduhaev and S. A. Kucherov, Phys. Lett. 112A, 61 (1985). 

111 A. D. Migone, Z. R. Li, and M. H. W. Chan, Phys. Rev. Lett. 53, 810 (1984); A. J. Jin, M. R. 
Bjustrom, and M. H. W. Chan, ibid. 62, 1372 (1989). 

112 k. J. Strandburg, Phys. Rev. B 46, 11190 (1992). 

113 88, 1622 (1985). 

11^ M. P. Ahen, D. Frenkel, and W. Gignac, J. Chem. Phys. 78, 4206 (1983); A. D. Novaco and 
P. A. Shea, Phys. Rev. B 26, 284 (1982); R. K. Kaha and P. Vashishta, J. Phys. C 14, L643 
(1981). 

115 V. N. Ryzhov, E.E.Tareyeva, Phys. Rev. B 51 N 14, 8789-8794 (1995). 

116 V. N. Ryzhov, E.E.Tareyeva, JETP (Moscou), 108, 2044 (1995). 

11'^ V.N.Ryzhov, E.E.Tareyeva, Theor. Math.Phys. (Moscou), 73, 463 (1987). 

118 V.N.Ryzhov, E.E.Tareyeva, J. Phys. C, 21, 819 (1988). 

119 V.N.Ryzhov, Theor. Math.Phys. (Moscou), 80, 107 (1989). 

120 V.N.Ryzhov, J. Phys.: Condens. Matter, 2, 5855 (1990). 

37 



V.N.Ryzhov, E.E.Tareyeva, Phys.Lett.A 158, 321 (1991). 

V. N. Ryzhov, Zh. Eksp. Teor. Fiz. 100, 1627 (1991) [Sov. Phys. JETP 73, 899 (1991)]. 
V.N.Ryzhov, E.E.Tareyeva, Theor.Math.Phys.(Moscou), 92, 331 (1992). 
V.N.Ryzhov, E.E.Tareyeva, Theor.Math.Phys.(Moscou), 96, 425 (1993). 

125 V. N. Ryzhov, E. E. Tareyeva, Physica C, 205, 55 (1993). 

126 A. B. Dzyubenko and Yu. E. Lozovik, Zh. Eksp. Teor. Fiz. 102, 284 (1992). 

127 b. J. Alder and T. E. Wainwright, Phys. Rev. 127, 359 (1962). 

128 J. A. Barker and D. Henderson, Rev. Mod. Phys. 48, 587 (1976). 

12^ J. Pearl, in Low Temperature Physics LT9, edited by J. G . Daunt, D. O. Edwards, F. J. Milford, 
and M. Yagub (Plenum, New York, 1965), p. 566. 

130 V.N.Ryzhov, E.E.Tareyeva, Phys.Rev.B, 1993, 48, 12907. 

131 V.N.Ryzhov, E.E.Tareyeva, Phys.Rev.B, 1994, 49, 6162. 
i32D.Yu.Irz, V.N.Ryzhov, E.E.Tareyeva, Phys.Lett.A, v.207, 374 (1995). 

i33D.Yu.Irz, V.N.Ryzhov, E.E.Tareyeva, Theor.Math.Phys. (Moscow) 104, 337 (1995); ibid 107, 
100 (1996). 

i34D.Yu.Irz, V.N.Ryzhov, E.E.Tareyeva, Phys.Rev.B. 54, 3051 (1996). 
135 J. M. P. Hsu and A. Kapitulnik, Phys. Rev. B 45, 4819 (1992). 



38 



